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Within the framework of proton model with taking into account the piezoelectric interaction with the shear 
strain eg, a dynamic dielectric response of KD2PO4 type ferroelectrics is considered. An expression for the 
piezoelectric resonance frequencies of the rectangular thin plate of the crystal cut in the (001) plane is 
obtained. 
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1. Introduction 

In our previous paper [ [T] we explored the dynamic dielectric response of square thin plates 
cut from the KH2PO4 family ferroelectric crystals in the planes (001), perpendicular to the axis 
of spontaneous polarization. Using the proposed modification of the proton ordering model [ [2] 
that includes the piezoelectric coupling with the shear strain eg, within the framework of the 
Glauber approach [ [3] and the four-particle cluster approximation, we obtained an expression for 
the dynamic dielectric permittivity of a crystal, which that takes into account the dynamics of the 
shear strain eg- In the low- frequency limit this expression coincided with the static permittivity 
of a mechanically free crystal, whereas in the microwave region it coincided with the dynamic 
permittivity of a clamped crystal, exhibiting a relaxational dispersion. 

In the intermediate region, the obtained permittivity had numerous peaks associated with the 
piezoelectric resonances. It should, however, be mentioned that while solving the partial differential 
equations for the strain in [[1], the boundary conditions were not set correctly. Instead of demanding 
that the entire edges of the plate were mechanically free, we considered the plate free only at its 
vertices. It resulted in the underestimated values of the resonant frequencies. In the present paper 
we shall correct this error. 

We shall not repeat here the details of the previous calculations, which were correct. The system 
Hamiltonian, most of the used notations, as well as derivation of the dynamic dielectric permittivity 
of a clamped crystal (the pseudospin subsystem dynamics), can be found in [[1]. 



2. Calculations 

We shall consider vibrations of a thin x Ly rectangular plate of a KD2PO4 crystal, cut in 
the (001) plane, induced by time-dependent electric field E^t = £"36*"*. In the ferroelectric phase 
this field, in addition to the shear strain Eq, induces also the diagonal components of the strain 
tensor Si, but for the sake of simplicity we shall neglect them. 

Dynamics of pseudospin subsystem will be considered in the spirit of the stochastic Glauber 
model [13], using the four-particle cluster approximation. The system of equations for the time- 
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dependent deuteron (pseudospin) distribution functions is 
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where £q-fi{t) is the local field acting on the f'th deuteron in the gth cell, which can be found 
from the system Hamiltonian (see [[T]); a is the parameter setting the time scale of the dynamic 
processes in the pseudospin subsystem. 

Taking into account the symmetry of the distribution functions 
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from we obtain a closed system of equations for the time-dependent single-particle, pair, and 
three-particle deuteron distribution functions 
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The used here notations can be found in [\T\. 

Dynamics of the deformational processes in KD2PO4 is described using classical Newtonian 
equations of motion of an elementary volume 
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where p is the crystal density, Ui are the displacements of the elementary volume along the axis 
Xi, is the mechanical stress. The shear strain eg is determined by the displacements ui and U2 
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The expression for the shear stress axy — erg, being the function of jy^^''^, and eg, is found from 
the constitutive equations derived in [[T]. 

At small deviations from the equilibrium we can separate in the systems ([3]) and ^ the static 
and time-dependent parts, presenting the dynamic variables as sums of the equilibrium values and 
of their fluctuational deviations 
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The fluctuational parts are assumed to be in the form of harmonic waves 

=,7«(x,y)e-S 4'^ ^ rj^'\x,y)e^-\ . . . 
uu = ui(2;)e'"S U2t = U2{x)e''^\ = ee{x, y)e^"*. 

The fluctuational part of (j3]) is reduced to the system of linear first-order differential equations 
with constant coefficients, solving which we get 
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the notations introduced here can be found in [[T]. 
Equations ^ reduce to 



■ klm = 0, ^ + klu2 = 0, 



where fcg is the wavenumber 



and 
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Differentiating the first and second equations of ((T]) with respect to y and x, correspondingly, 
remembering that we neglect the diagonal strains ei — dui/dx and £2 = dui/dy, and adding the 
two obtained equations, we arrive at the single equation for the strain eg 



dy^ 



+ klee{x,y) = 0. 
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Boundary conditions for £6(2:, y) follow from the assumption that the crystal is simply supported, 
that is, it is traction free at its edges (at x ~ 0, x ~ L^, y — Q. y ~ Ly, to be denoted as E) 



o-els = 0. 
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In our previous consideration [ Ij this condition was fulfilled at the corners of the crystal plate 
only, not at all its edges. Substituting pT|) into the constitutive relations, we obtain the explicit 
boundary conditions for the strains in the following form 
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where 
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Solution of ([8]) with the boundary conditions ([12]) is 



£G{x,y) = £60 + £60 X! 
with Lo^i given by 
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Using the expression, relating polarization P3 to the order parameter t^*-^) and strain £g (see [ 
[T]), we find that 

P^(x,y,t) = P^{x,y)e^\ (16) 
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where 
and 

is the dynamic dielectric susceptibihty of a clamped crystal. 

Now we can calculate the dynamic dielectric susceptibility of a free crystal ^33(0!'^) 
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From (fTH. we find that 
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Finally the dynamic dielectric susceptibility is 

x?i(c.) = x?i(^) + i?6(^)4r^- (19) 

In the static limit (w — >■ 0, i?6('^) ^ 1) from (fTO|) we obtain the static susceptibility of a free 
crystal [%■, in the high frequency limit (Em=o 64/[(2A: + 1)^(2^ + ifir'^] = 1, and Re{uj) 0) we 
get a dynamic susceptibility of a mechanically clamped crystal, exhibiting relaxational dispersion 
in the microwave region. Thus, eq. (|19p explicitly describes the effect of crystal clamping by high- 
frequency electric field. 

In the intermediate frequency region, the susceptibility has a resonance dispersion with numer- 
ous peaks ar frequencies where Re[i?6('^)] ^ oo- Frequency variation of c^g(aa;) is perceptible only 
in the region of the microwave dispersion of the dielectric susceptibility. Below this region c^g(aa;) 
is practically frequency independent and coincides with the static elastic constant cf^. Since the 
resonance frequencies are expected to be in the 10'* — 10^ Hz range, depending on temperature 
and sample dimensions, we can neglect the frequency dependence of c^g(Q!aj) ^ and reduce the 
equation for the resonance frequencies (1151) to an explicit expression by putting in it cfg{auj) — > c^g. 

Comparing psp to the expression obtained previously [ 1 for a square L x L plate cut in the 
(001) plane ' 

7r(2fc + l) /cfg(cjfc) 

we can see that the incorrectly set boundary conditions [ 1] led to the times smaller lowest 
resonance frequency than the correct one. However, the low and high frequency limits of the 
susceptibility (the static value and the clamped values with the relaxational dispersion in the 
microwave region) were correct. 

3. Conclusions 

Within the proton ordering model with taking into account the shear strain eg we explored a 
dynamic response of the KD2PO4 type crystals to an external harmonic electric field E3. Dynamics 
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of the pseudospin subsystem is described within the stochastic Glauber approach. Dynamics of the 
strain eg is obtained from the Newtonian equations of motion of an elementary volume, with 
taking into account the relations between the order parameter of the pseudospin subsystem and 
the strain in the static limit. A corrected expression for the piezoelectric resonance frequencies of 
the rectangular thin plates of these crystals cut in the (001) plane is obtained. 

The ultimate goal of the present studies is will be to generalize the obtained expression for 
the dynamic permittivity the case of the mixed Rbi_2;(NH4)a;P04 type systems, exhibiting the 
proton glass properties, in order to explore the dynamic dielectric response. It is known that, just 
like their pure constituents, these mixed systems are piezoelectric, and their dynamic dielectric 
permittivity has a piezoelectric resonance dispersion [|4]. As our preliminary calculations show, 
the experimentally obtained resonant frequencies of such mixed crystals [ 5] are well described by 
the obtained here expression for the resonant frequencies, provided that the corresponding elastic 
constant cfg of such a system is known. 
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